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Abstract 

The Gaudin model based on the sZ2-invariant r-matrix with an extra 
Jordanian term depending on the spectral parameters is considered. The 
appropriate creation operators defining the Bethe states of the system are 
constructed through a recurrence relation. The commutation relations 
between the generating function t(X) of the integrals of motion and the 
creation operators are calculated and therefore the algebraic Bethe Ansatz 
is fully implemented. The energy spectrum as well as the corresponding 
Bethe equations of the system coincide with the ones of the sZ2-invariant 
Gaudin model. As opposed to the sZ2-invariant case, the operator f (A) and 
the Gaudin Hamiltonians are not Hermitian. Finally, the inner products 
and norms of the Bethe states are studied. 
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1 INTRODUCTION 



1 Introduction 

A model of interacting spins in a chain was first considered by Gaudin [TJ [2] . 
In his approach these models were introduced as a quasi-classical limit of the 
integrable quantum chains. Moreover, the Gaudin models were extended to any 
simple Lie algebra, with arbitrary irreducible representation at each site of the 
chain [2]. 

In the framework of the Quantum Inverse Scattering Method (QISM) [3]- 
[5] the Gaudin models are based on classical r-matrices which have a unitarity 
property 

r\2(n,v) = ™?"2i0,m)- (1-1) 
Then the Gaudin Hamiltonians are related to classical r-matrices 

N 

=J2 r «b(z a ,z b ). (1.2) 

bjia 

The condition of their commutativity [H^ a ',H^ b '] — is nothing else but the 
classical Yang-Baxter equation 

[r a b{z ai z b ), r ac (z a ,z c ) + r bc (z b , z c )] + [r ac (z a , z c ), r bc (z b , z c )] = (1.3) 

where r belongs to the tensor product g eg) g of a Lie algebra g, or its repre- 
sentations, and the indices fix the corresponding factors in the TV-fold tensor 
product. 

Considerable attention has been devoted to Gaudin models based on the 
classical r-matrices of simple Lie algebras [S]-[II] and Lie superalgebras [T51 [TBI 
IT7] . The spectrum and corresponding eigenfunctions were found using different 
methods such as coordinate and algebraic Bethe ansatz, separated variables, 
etc. Correlation functions for Gaudin models were explicitly calculated as com- 
binatorial expressions obtained from the Bethe ansatz [3] . In the case of the s?2 
Gaudin system, its relation to Knizhnik-Zamolodchikov equation of conformal 
field theory [14] or the method of Gauss factorization [TTJ, provided alterna- 
tive approaches to computation of correlations. The non-unitary r-matrices and 
the corresponding Gaudin models have been studied recently, see [IB] and the 
references therein. 

A classification of solutions to classical Yang-Baxter equation fj 1 . 3[) . with 
the property (jl.ip . has been done in [19]. In particular, a class of rational 
solutions corresponding to the simple Lie algebras sl n was discussed in [19] and 
then further studied and generalized in [501 IS]- In the s h case, up to gauge 
equivalence, only two such solutions exist: the s^-invariant r-matrix with an 
extra Jordanian term (deformation) 

r® 

r({^v) = -^-+^(h®X+ - X+ ®h), (1.4) 

where cf = h ® h + 2 {X + ® X~ + X~ ® X + ) denotes the quadratic tensor 
Casimir of s/2, and the s?2-mvariant r-matrix with a deformation depending on 
the spectral parameters 

c® 

r e (/*, v) = —2— + £ (h O (vX+) - (nX+) ® h) . (1.5) 
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The Gaudin model based on the r-matrix (|1.4p was studied as the quasi- 
classical limit of the corresponding quantum spin chain |221 123] . which is a 
deformation of the XXX spin chain where the Jordanian twist is applied to the 
Yang R-matrix [MJ [35] . The algebraic Bethe Ansatz for this Gaudin model was 
fully implemented in [26] , following ideas given in (TSJ US [27] . 

In the framework of the QISM, the classical Yang-Baxter equation guaran- 
tees for the Gaudin model the existence of integrals of motion in involution. A 
more difficult problem is to determine the spectra and the corresponding eigen- 
vectors. The Algebraic Bethe Ansatz (ABA) is a powerful tool which yields 
the spectrum and corresponding eigenstates in the applications of the QISM 
to the models for which highest weight type representations are relevant, like 
for example Gaudin models, quantum spin systems, etc [51 [5J [TO] . However the 
ABA has to be applied on the case by case basis. 

An initial step in the application of the Algebraic Bethe Ansatz to the 
Gaudin model based on the r-matrix (jl.5D . is to define appropriate creation 
operators which yield the Bethe states and the spectrum of the corresponding 
Gaudin Hamiltonians. The operators Z?j^ (//i, ■ ■ ■ ;A*m) used for this model, in- 
troduced in [28) . are defined through a recursive relation that generalizes the 
one used in [2 6) for the Gaudin model based on the r-matrix (jl.4p . The com- 
mutation relations of the ■ . ■ , Hm) operators with the entries of the 
L-operator are straightforward to calculate, see lemma 12.21 The main step, 
lemma 12.31 below, in the implementation of the ABA is to calculate explicitly 
the commutation relations between the creation operators Bm(^i, ■ • ■ , A*m) and 
the generating function of integrals of motion t(X) of this Gaudin model. Fi- 
nally, it is confirmed that the operator t(X) is not Hermitian with respect to the 
natural inner product. 

Besides determining the spectrum of the system, via ABA, another impor- 
tant problem is to calculate the inner product of the Bethe vectors. To this end 
it is necessary to introduce the duals of the creation operators _B^j(/ii, . . . , (j,m)< 
A method developed in 26 to determine the inner product of the Bethe vectors 
will be used here. The main idea is to establish a relation between the creation 
operators of the system and the ones in the sZ2-invariant case, i.e. when £ = 0. 
The advantage of this approach is that the cumbersome calculation involving 
strings of commutators between entries of the L-operator becomes unnecessary. 
As opposed to the sZ2-invariant case, the Bethe vectors ^ m (p-i, ■ • • , Mm) for the 
deformed Gaudin model are not orthogonal for different M's, with respect to 
the natural inner product considered. 

The article is organized as follows. In section 2 we study the deformed 
Gaudin model based on the r-matrix (ll.5[) , emphasizing the operators B M and 
their properties. Using the commutation relations of the creation operators Bm 
with the generating function of integrals of motion t(X), in section 3 the Bethe 
vectors of the system are given. It is also confirmed that the spectrum is the 
one of the invariant system. The duals of the creation operators are introduced 
in section 4 in order to obtain the expressions for the inner products and the 
norms of Bethe states. The results of the study are briefly summarized in the 
conclusions. The proofs of the lemmas are given in the appendix. 
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2 Algebraic Bethe Ansatz for deformed Gaudin 
model 

The Gaudin model based on the s^2-hivariant r-matrix with an extra Jordanian 
term depending on the spectral parameters 



r 5 (A,/i) 



A — /j, 



1 



h®h + 2(X^ 



A — /j, 

+ nh®(nx+) - (ax+) 



X 



x- 



,*+) 



where are the standard s?2 generators 

[h,X + ] = 2X+, [X+,X-]=h, 



[h,x- 



-2X~ 



(2.1) 



(2.2) 



is considered. The matrix form of rj(A, fi) in the fundamental representation of 
s?2 is given explicitly by [35] 



/- 



1 



r s (A,/x) 



A — (i 




V 





-A£ 


1 


2 


A — /i 


A — /i 


2 


1 


A — /! 


A — yU 









\ 

A£ 

1 



(2.3) 



here A, /i G C arc the so-called spectral parameters and £ G C is a deformation 
parameter. The next step is to introduce the L-operator of the Gaudin model 
[9]. In the present case the entries of the L-operator 



L(X) = 



h(X) 
2X+(X) 



2X-(xy 

-h(X) i 



(2.4) 



are given by 



N 



MA) = J2 



N 



a=l 



X: 



X - z a 



X — z, 

-Xh n 



■iz a xt 



N 



a=l 



X — z. 



(2.5) 



MJ 



where h a = n^\h) G End(K ( a> ), X± = ^i M (X±) G End(V^ J ) and 4' 
is an irreducible representation of s/2 whose representation space is Va a cor- 
responding to the highest weight t a and the highest weight vector oj a G V a , 
i.e. X+ui a = and h a ui a — £ a u a , at each site a = 1, . . . , N. Notice that £ a is 
a nonnegative integer and the (£ a + l)-dimensional representation space Va 
has the natural Hermitian inner product such that (cf. |29| ) 



x: 



(X-y=X+ and h* a 
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The space of states of the system % — V± (g> ■ ■ ■ ® V^' N%> is naturally equipped 
with the Hermitian inner product (•(•) as a tensor product of the spaces VJj 
for a= 1,...,N. 

The L-operator (|2.4[) satisfies the so-called Sklyanin linear bracket [9] 



r^X,fi),L(X)+L(fj,) 



(2.6) 



Both sides of this relation have the usual commutators of the 4x4 matrices 
L(A) = L(A)®1, L(X) = 1(E) L(X) and r^(X,/j,), where 1 is the 2 x 2 identity ma- 
trix. The relation (|2.6[) is a compact matrix form of the following commutation 
relations 

[h(\),hQi)]=2£ (XX + (X)-fxX+(fi)) 

[x-(x),x-(»)] = -e (»x-(x) - \x- (ji)) , 



[x+(x),x-(n)]=- 

[X+(\),X+M] =0, 
[ft(A),X-(^)] 



fe(A) - h(n) 
X — /J, 



+ ^X+(X), 



(2.7) 



X-(X) -X-(u) . i / \ 
2 — — ^+£MG«), 



^(A),X + (/i)] =-2 



X — n 
X+(X)-X+^) 
X — n 



In order to define a dynamical system besides the algebra of observables a 
Hamiltonian should be specified. Due to the Sklyanin linear bracket (|2.6p the 
generating function [9] 



t{X) = - TrL 2 (A) = h 2 (X) - 2h'(X) + 2 (2X~(X) + £A) X + (X) (2.8) 
satisfies 

t(X)t(p)=t(jj,)t(X). (2.9) 

The pole expansion of the generating function t(A) is obtained by substituting 
HHg) into d^SJ 



A' 



*w = E 



es=l 



2i/( Q ) 



(A - z a ) 2 A - z a 



N 



.,6=1 



(2.10) 

The residues of the generating function t(X) at the points A = z a , a = 1, . . . , N 
are the Gaudin Hamiltonians 

H[a) = £ /c^6) + ^ (2b ^ x + _ ZaX +^ fc) \ ; (2 n) 

where 02(0, b) = h a h\, + 2(X+X^ + X~X^). In the constant term of the 
expansion (12.101) the notation 



N 



Y, 



(<//) 



5>, 

a=l 



(2.12) 
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for Y = (hfX^), was used to denote the generators of the so-called global sZ 2 
algebra. In the case when £ = the global s/2 algebra is a symmetry of the 
system. 

Moreover, from the formula (I2.10[) it is straightforward to obtain 

t(A) - t(A)„ + 2£ (h(X) X+ gl) + X+ gl) - Xh (gl) X+(X)) + e(X+ gl) ) 2 , (2.13) 

where X^ gV) = Y,a=x z aX+ , h(X) = h(X)\ i=0 and t(X) = t(X)\ i=0 is the gener- 
ating function of the integrals of motion in the s^-invariant case. 

In order to obtain the Bethe states of the system it is necessary to fully 
implement the ABA. In the space of states T-L the vector 

tt + = lui <8> • • • <g> uj n (2-14) 

is such that (f2+|f2+) = 1 and 

A+(A)0+ = 0, h(X)ft + = p(A)Q+, (2.15) 

with 

P(A)=EA- (2-16) 

The commutators between the generating function (|2.8[) and the entries of 
the L-operator (|2.5p are 

[t(A),JC+(M)] =4 W)* + (M)-MM)* + (A) - H XX+(X)X+(,), (2.17) 
. 8 X-MX^(A)-A-(A)X + (,.) _ 4?(1 + , MA))X+W , (2 , 8) 

A — /I 

mx-M] =4 ^-(A)MM)-x- W MA) +2g(1+M(A))%) ^ 

+ 4eAA-(/i)A+(A) + 2(^) 2 A+( M ). 

The next step in the Algebraic Bethe Ansatz is to define appropriate creation 
operators that yield the Bethe states. The creation operators used in the SI2- 
invariant Gaudin model coincide with one of the L- matrix entry [21 19]. However, 
in the present case these operators are non-homogeneous polynomials of the 
operator A - (A). It is convenient to define a more general set of operators in 
order to simplify the presentation. 

Definition 2.1. Given integers M and k > 0, let fi = {fii, ■ ■ ■ , /jm} be a set of 
complex scalars. Define the operators 

BSV) = (*~(Mi) + (X-(m) + (k + 1)^2) ■ • ■ (A-(/i M ) + (M + k - 1K M m) 

M+k-l 

= Yl (X-(n n - k +i)+n^ n - k+1 ) , (2.20) 



with B ( k) = 1 and B$ = for M < 0. 
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The following lemma describes some properties of the B^(fi) operators 
which will be used later on. 

Lemma 2.1. Some useful properties of B^il^i, . . . ,/Um) operators are: 

i) ^mVi, • • - ,A»m) = B${m, . . . ,/x m )B^+™ ) (^ m+ i,. . . ,/i M ), m = 1,...,M 
= i?f ) ( Ml )<+ 1 1 ) ( M2 ,..., MM ) 
= ^m-i^i. • • • ,MM-i)M M+fe_1) (MM); 

(2.21) 

ii) TTie operators B^(/j,i, . . . ,/xm) are symmetric functions of their argu- 
ments; 

hi) b2V) = sr i) (M)+cE^rf_ >«). 

The notation used above is the following. Let /x = . . . , ^m} be a set of 
complex scalars, then 

/i (<1, - ,ifc) =/x\{/x il ,..., W J 

for any distinct . . ,ik £ {1, . . . , M}. 

The Algebraic Bethe Ansatz requires the commutation relations between the 
entries of the L-operator and the B^ (fi) operators. To this end the notation 
for the Bethe operators are introduced 

M (\-p) = h(\) + Y,-zrx> ( 2 - 22 ) 

where fi = . . . , hm-i} and A € C \ /z. In the particular case M = 1, 
/3i(A;0) = /i(A) and will be denoted by /3i(A). The required commutators are 
given in the following lemma. 

Lemma 2.2. T/ie commutation relations between the operators h{\), X ± (X) 
and the B^{n\, . . . ,/jm) operators are given by 

<=i A ~ Mi 

+ ^E B m-i ) ^ W )( v ^(^;^ ) ) -2fcj; (2.23) 



M R (fe+l)/i, I | „(»,j)l 



i<j 



M 



/3 M (A;/xW)-/3 M (^;/xW) 



A - Mi 

(2.24) 
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M 
i=l 

M 

= B^ +1 (X U M ) - e E ^ B m ( A u (2-25) 

i=l 

It is important to notice that the creation operators that yield the Bethe 
states of the system are the operators B^ (fi), below denoted by £?m (/•*)■ These 
so-called B-operators are essentially based on those proposed by Kulish in [27] 
and then used in the implementation of the ABA for the Gaudin model based 
on the r-matrix (jl.4[) [2B]. A recursive relation defining the creation operators 
follows from i) of the lemma |2~T1 

B M (p) = Bm-i(m (M) ) (X-(h m ) + (M - l)£ii M ) . (2.26) 

As a particular case of lemma [2~2l the commutation relations between operators 
h(X), ^ ± (A) and the B-operators can be obtained by setting k = in (12. 23|) - 
(12351) . 

The main step in the implementation of the ABA is to calculate the com- 
mutation relations between the creation operators, in this case the B-operators, 
and the generating function of the integrals of motion i(A). The complete ex- 
pression of the required commutator is given in the lemma below and is based 
on the results established previously in this section. This will be crucial for 
determining the spectrum of the system and thus constitutes one of the main 
results of the paper. 

Lemma 2.3. The commutation relations between the generating function of the 
integrals of motion t(X) and the B-operators are given by 



t(X)B M (n) = B M (n) lt(X) - J2 + E n Ul V + 4M£AX+(A) 

A ^ Bm(XU^) 5 . m , 

M 

+ 2^E B M-i(M (4) )(M i MA) + i)$ m (jm;^) 

1=1 

+ 4C 2^ w r~ pm{ih;h k> ) 

M 
M 



2e 2 E^M-i(M W )^ + (^)- (2-27) 



Proof. In the case when M = 1 the commutator between the operators t(X) and 
Bi(n) = X~(jjl) has been calculated already, see (|2.19p . 
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The next step is to outline the proof in the case when M > 1. The starting 
point in the calculation of the commutator [t(X), Bm (m)] has to be the expression 
so that 



[t(\),B M (n)) = [/.(A), [h(X),B M (fi)]] + 2 [h(X),B M (jj)] h(X) - 2^ [h(X),B M (p)] 

+ 2(2X-(A)+£A) [X+(A),S M (/x)] + 4 [X~(X), B m (m)] * + (A). 

(2.28) 

Notice that the terms on the right hand side of the equation above involve only 
the commutators between the entries of the L-operator and the B-operators. 
Each of them will be calculated separately, but previously an auxiliary result is 
established. Then the first term on the right hand side of the equation (|2.28[) 
can be calculated using the expression (|2.23[) 

M 

[h(X), [/1(A), S M (/X)]] = B M (jj) ]T J- r~ r 



^ /Bm(AUm«)-BmM ^B A/ (AU/xW) N 



(A-/i,) 2 A-/Xj 



2=1 \ 

^ (A - m») (A - Mj) ^ A -/Xi ^ A -^ 
+ 4 ^E(A-, J )(A-M J ) +2 %^- l(/X } A^ 

+e|> (4 ^- l(Au ^ ) :^- l( ^ )) +^ 2 (^) (m./Wm;;^) - 2)) 

Af 

+ 2£ 2 ^>5«_ 1 (m«)(AX+(A) - M i* + (Mi))- (2.29) 

i=i 

The equation (I2.23[) . with k = 0, is also used to determine the next two terms 
on the right hand side of the equation (12.28)) 

2 [MA), B„W\ h(X) = B„W J2 + 4 £ ^ ( l A}U/iW) MA) 

M M 
i=l i=l 

(2.30) 

and after differentiating (I2.23[) . with k = 0, 

- 2 A [M A), B„ („)] - 4 1 ^'^r'"' B """ - AJ; " J(AUM '" ) - 
dA ^ (A - Mi) A - Mi 

(2.31) 

The forth term on the right hand side of the equation (12.28)) can be obtained di- 
rectly using the expression for the commutator [X + (X), Bm(m)] given in (|2.24l) . 
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with k = 0, 

2(2X-(A)+£A) [X+(A),%(M)] = -4^%(AU^)^^M^ 
g ^ X-(A) B^jX U /x^)) ,y g( i) , (i) J M (X-,pW)-& M (iHwW) 

f£ (\-m)(x-Hj) M_1 a-/h 

M 

+2e\Y J ^ B M-M [t) ) x+ w- ( 2 - 32 ) 

Finally, the last term on the right hand side of the equation (|2.28p is determined 
using the equation (|2.25[) . with k = 0, 

M 

4[X-(X),B M (n)] X+(X) =4M^ASM(/x)X+(A)-4^ Mi S M (AU/LX»)X+(A). 

i=l 

(2.33) 

The last step in the proof of the equation (|2.27[) is a straightforward addition 
of all the necessary terms (|2.29[) - (|2.33[) calculated above. □ 

Computationally the equation (|2.27[) is the crucial step of the ABA. In the 
next section this result will be used to determine the spectrum of the system. 

At the very end of this section a property of the B-operators relevant for the 
calculation of the inner products and the norms of the Bethe states is estab- 
lished. Namely, using the realization (|2.5[) . the relation between the B-operators, 
given by the recurrent relation (|2.26p . and the corresponding ones for the sfo- 
invariant Gaudin model, the case when £ = 0, is obtained 

Af-l M 



B M (n) =Y,t k E B M _ k (» jl ,...^ jM _ k ) pl M -V (^'•••»^-* ) )+£ M i5$(/*) 

(2.34) 



k = Q h<-<3 M -k 

where 



and the operators p^(jii, . . . ,/Xj) are given by 

p»( m ,..., Mi )=Mi(-^+J/')---M i (-^+i+i-l), and p« = 1. (2.35) 

In the next section it will be shown how the full implementation of the 
Algebraic Bethe Ansatz is based on the properties of the creation operators 
Bm(/i) obtained here. 



3 Spectrum and Bethe vectors of the model 

In this section the spectrum and the Bethe vectors of the Gaudin model based on 
the r-matrix (|1.5[) will be obtained by the ABA. The first step of the ABA is the 
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observation that the the highest spin vector f2+ (|2.14l) is an eigenvector of the 
generating function of integrals of motion t(X), due to (|2.8[) . (|2.15p and (|2.16p . 
Then Bethe vectors ^m(aO are given by the action of the -B-operators on the 
highest spin vector \1/m(m) = Bm{h)Q+- Finally the spectrum of the system 
is obtained as a consequence of the commutation relations between t{X) and 
Bm(h) lemma [2~31 (I2.27[) . The unwanted terms coming from the commutator 
are annihilated by the Bethe equations on the parameters fi — {pi, . . . , pm} 
as well as the condition X + (X)il + = 0. Thus the implementation of Algebraic 
Bethe Ansatz in this case can be resumed in the following theorem. 

Theorem 3.1. The highest spin vector f2 + (|2.14p is an eigenvector of the op- 
erator t(A) 

t(x)Q + = A (\)n + (3.i) 

with the corresponding eigenvalue 

MA, . AA) - W) - £ jh: (t ^t) + 1 ■ M 

a=l \b=£a I o=l v ' 

Furthermore, the action of the B-operators on the highest spin vector VL + yields 
the Bethe vectors 

* M {n) = B M {iJL)n + , (3.3) 

so that 

t(X)^ M (p) = A M (A;/x)* M (/x), (3.4) 

with the eigenvalues 

dp M ° 

A m (A;m) = Pm(A;m) - 2— — (A;/i) and p M (A; /i) = p(A) - ^ 



dA r— X — Hi 

(3.5) 

provided that the Bethe equations are imposed on the parameters fi — {fJ,±, . . . , /in} 

pm(m»;^ (i) )-E^^-E^— = ' ' 1 ^ (3-6) 

o=l Ml Zq i^i Ml Mj 

Proof. The first step in the proof of the Theorem l3.1l is to show that the highest 
spin vector fi_|_ (|2.14p is an eigenvector of the generating function of integrals of 
motion t(X). Using ([2T5]> . (12TT5)) and (12~5|) one gets 

i(A)0+ = (/i 2 (A) - 2/i'(A) + 2 (2X-(X) + £A) X+(X)) n + = A (A)ft + 

with the eigenvalue Ao(A) = (p 2 (A) — 2p'(A)) and the function p(X) is given in 
the equation (|2.16p . The next step is an important observation that the action 
of the operator i(A) on the Bethe vectors ^m(/-*) (|3.3p is given by 

i(A)# M (M) = t(X)B M (fj,)n + = A (A)*m(m) + [*(A), B M (p)] «+• 

From the equation above it is evident that the main step in the calculation of 
the action of the generating function of the integrals of motion on the Bethe 
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vectors is computing the commutator between the operator t(A) and the creation 
operators Bm{h)- From the Lemma \2 . 31 follows 

t(A)<M/x)= [A a (X)-J2P^ + J2n Wl 

I f=t A ~^ ^ (A-Mij(A-Mi) / 

+ 4^ M M ) pm(m,;m w ) + 2^^(^;^ ) )x 



A — /i 

M 

x[ (l + Mi p(A))< ) _ 1 ( / iW) + 2 Ati ^; 

M / M 



(3.7) 



The unwanted terms in (|3.7p are all proportional to the function pu^i\ fJ-^) 
and thus vanish when the Bethe equations Pm (/•*»> M ) = 0' <|3.6p are imposed 
on the parameters /x = {/Xi, . . . , /jm}- Then the equation (|3.7p has the required 
form (|3.4p with the eigenvalues given by formula (|3.5p . □ 

Although the deformed Gaudin model does not admit the global SI2 sym- 
metry, it is of interest to notice that, in the cases when M could be equal to 
1 + ^^2a=i^a, the on-shell Bethe vectors (I3.3P are highest weight vectors for 
the action of the global SI2 algebra on H i.e., 



1 N 

X+ gl) ^ M (fi) = and h {gl) y M (p) = (-2M)* M (/i), when M = 1 + - ^4 



0=1 



and the Bethe equations (I3.6p are imposed on the parameters /x = {/xi, • ■ • , Mm}- 
In particular, the formulae (|2. 51) yield 



X M) = , lim ( A * + ( A )) and ^) = ( -7) , lim f^T^ ) • ( 3 - 8 ) 



L W ~ A-too ^ " w ~ V U x-+°a V A 

The statement above follows from the formulae (I3.8P , the lemma 12.21 the equa- 
tions (j2~2lT) and (j2~2"5j) . and the lemma I2TT1 the equation (|2~2Tj) . 
From the commutator 

+ 2^ (s0 + e (h w + h {gl) ) X+ gl) + 2ex+ iy (3.9) 

it follows that, even if the dimension of the invariant subspace associated with 
each Bethe vector is the same as for the invariant model, the vectors (X^ gl ^) n ^ m(/-0 
are not eigenvectors of i(A) and, in general, they even do not generate the cor- 
responding invariant subspace. 

To obtain the spectrum of the Gaudin Hamiltonians H ^ (|2.1ip one uses the 
pole expansion (|2.10p of the generating function i(A) and the previous theorem. 
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3 SPECTRUM AND BETHE VECTORS OF THE MODEL 



Corollary 3.1. The Bethe vectors ^/(/-O (|3-3[) are eigenvectors of the Gaudin 
Hamiltonians (|2.11[) 

H^^ M (fj,) = E$* M (fi), (3.10) 
with the corresponding eigenvalues 

N o o M Of 

b=£a i=l 

Below it is assumed that the local parameters z a are real numbers, i.e. 
z a = z a for all a — 1, . . . , N and the Hilbert structure in the representation 
space T-L is as in section 2. As expected, the spectra of the Gaudin Hamiltonians 
H ( a ) are real, although these operators are not Hermitian, for all a = 1, . . . , N, 

2f(«0 / (H^T = E ( + ZfahaX; - ZahbX^)) . (3.12) 

\ z a z b J 
b^a x ' 

Analogously, the generating functions of the integrals of motion t(X) is also not 
Hermitian 

t(X) * (t(A))* = t(X) +2C (x w h(X) Q + X w -\X-(X)h {gl) - ^hl gt )+e(X w )\ 



where X ( l) = J2a=i z * x ~a 



(3.13) 



The problem of correlation functions of spin chains have been extensively 
studied [B] as well as for the sZ2-invariant Gaudin model [TT]. In the case of 
the inner product of the off-shell Bethe vectors, for the system under study, a 
similar approach to that in [26] is developed below. 

In order to establish a relation between the inner products of the Bethe 
vectors (|3.3[) and the corresponding ones in the s/2-invariant case it is of interest 
to consider the property (I2.34[) of the B-operators, established at the end of the 
section 2, in its dual form 

M-l M 

B* M W = E Z k E Pi M ~ k) B* M _ k (X il ,...,X iM _ k )o+i M p^ ) (X) 

(3.14) 



A.' — *!<-■-<; 



M - k 



where B* M {\) = (B M (X)) and 
B* N (M, ■ ■ • , AaOq - B* N (X U . . . , X N )\ i=0 = X+(X 1 )X+(X 2 ) ■ ■ ■ X+(X N ). (3.15) 
Notice that the operators (13. 14)) do not annihilate the highest spin vector 

B M (A 1; • • • , Am)«+ = S M f$(\)fl + = £ M PM (X 1 , . . . , X M )il + , (3.16) 
here p M (Xi, Am) = Ax • ■ • X M Pm with 

^=(-^f)(-^f + l)-(-^ + M-l) and „ W) =2> 

(3.17) 
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3 SPECTRUM AND BETHE VECTORS OF THE MODEL 



The following lemma yields an explicit relation between the inner prod- 
ucts of the Bethe vectors, when the Bethe equations are not imposed on their 
parameters ("off-shell Bethe states"), of the deformed Gaudin model and the 
corresponding ones in the s/2-invariant case. 

Lemma 3.1. Let Mi, M2 be two natural numbers and A = {Ai, . . . , AmJ, A* = 
{/ii, . . . , /ij\/ 2 } be sets of complex scalars. Then the inner product between the 
off-shell Bethe states (|3.3p is given by 

Mi-i A/2-1 

(*M 1 (A)|fM 2 (/x))=f Ml C M2 PM 1 (A)p M2 (/x)+ ]T f fcl £ fca ((p Ml , fcj ( A) 1 4>M 2 , fc 2 ) ) ) 

fc 1= fe 2 =0 



(3.18) 



where 

M 



(A) = £ A M - k) (A ' ' ) (Ai 1 ,...,Ai M _ fe ) , 

(3.19) 



l^il<---<iA/_ fc 

wzi/i *m(-)o = *m(-)I^=o awrf 

pp ) (A 1 ,...,A l )^ + =P?' ) (A 1 ,...,A l )r! + 

= (Ar(-f + i) • • • \ii- P -f + i + 3 ~ 1)) n+. (3 ' 20) 

Proof. Let ^E'mi(A) and V^m^m) be two off-shell Bethe states. Then, 

(V Ml (\)\V M2 (u)) = (S Ml (A)0+ I B M2 (n)n + ) = (0+ I B* Ml (\)B M2 (n)n + ). 

(3.21) 

From the eqns. (j2lHl) and (l3~TI)) it follows 

(X)B M2 (/x) = p<g (X)B M2 (n) + e h B* Ml (X)p<£ (/x) - f M ^ Ma (A)p£ 2 M 

A/i-1A/ 2 -1 Aifi A/ 2 



E E^ fc2 E E ^-^(^ 

fei=0 fe 2 =0 il<— <iMi-fci Oi<---<3M 2 ~k 2 



X 



x SM 1 - fcl (A ll ,...,A 4A/i _ fc JoBA/ 2 - fe2 (^ 1 ,-.-,^M 2 - fc2 )o4f 2 " fc2) (/' " ■ ) ■ 

(3.22) 

Substituing (|3.22[) into (|3.2ip . and taking into account (|3.16[) 

A/i-l A/ 2 -l 

(*A/i(A)|*M 2 (^)=p /l £ M2 PA/i(A)p M2 ( At )-|- ]T E X 

fc 1= k 2 =0 

A/i A/q 

X E E (A (il '-' iMl - fcl)) ) pif 2 - fe2) ^fai.-^-*,)) x 

il<---<iA/i-fci jl<---<jM 2 -k 2 

x /*M!-fei (A n , . . . , Ai Ml _ (5i ) *A/ 2 -fc 2 (Mil . • ■ • 7 HM 2 - k2 )o V 

(3.23) 
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Given that (-|-) is the Hermitian inner product in the representation space H, 
as considered in section 2, the right hand side of (|3.23[) can be expressed as in 
(13351) . □ 

The main result of this section is the equation Q3.18[) of the lemma 13.11 
where the inner products of the off-shell Bethe states of the deformed Gaudin 
model are expressed in terms of the s^-invariant inner products calculated by 
Sklyanin [TT]. A substantial advantage of this approach is that the result was 
obtained without doing the cumbersome calculation involving the commutators 
(|2.7[) between the entries of the L-operator. 

Given that (</>Mi,fci (A) |0M 2 ,fc 2 (/•*)) = unless M%—ki — M-i — k^-, the lemma 
13.11 yields an expression of the norm of the Bethe vectors. 

Corollary 3.2. The norms of the off-shell Bethe states ^ m{i^) are given by 

M-l 

\\*M(P)\\ 2 = \e M \v-?{Pm? + E I^UmA^W 2 (3-24) 

k=0 

where \n\ 2 = Mi| 2 ' ' ' I Mm 1 2 and </>m,/c(a*) are defined in p.!9[) . 

The corollary above reveals the full strength of our approach in the fact that 
the expression of the norm of the off-shell Bethe states again involves only the 
sZ2-invariant inner products given in |11] , 



4 Conclusions 

The Gaudin model based on the r-matrix (11.5[) is studied. The usual Gaudin 

realization of the model is introduced and the B-operators B^ (/ii, . . . , /im) 
are defined as non-homogeneous polynomials of the operator X~{\). These 
operators are symmetric functions of their arguments and they satisfy certain 
recursive relations with explicit dependency on the quasi- momenta fit, . . . , /jm- 
The creation operators Bm^i, ■ ■ ■ , (J-m) which yield the Bethe vectors form their 
proper subset. The commutator of the creation operators with the generating 
function of the Gaudin model under study is calculated explicitly. Based on the 
previous result the spectrum of the system is determined. It turns out that the 
spectrum of the system and the corresponding Bethe equations coincide with 
the ones of the s^-invariant model. However, contrary to the s/2-invariant case, 
the generating function of integrals of motion and the corresponding Gaudin 
Hamiltonians are not Hermitian. The explicit form of the generalized Bethe 
vectors associated to the Jordan canonical form of the generating function t(X) 
remains an open problem. 

In order to obtain expressions for the inner product of the Bethe vectors the 
operators B* m (ijl) are introduced. The method used in the study of the inner 
product of the Bethe vectors is based on the relation between the creation oper- 
ators of the system and the ones in the s/2-invariant case. The advantage of this 
approach is that the cumbersome calculation involving strings of commutators 
between the entries of the L-operator becomes unnecessary. The final expres- 
sion of the inner products and norms of the Bethe states is given in terms of 
the corresponding ones of the s/2-invariant Gaudin model. As opposed to the 
s/2-hivariant case, the Bethe vectors ^M(fi) for the deformed Gaudin model are 
not orthogonal for different M's. 
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Appendix 

A Proofs of lemmas 

The proofs of the lemmas, not included in the main text, are detailed in this 
appendix. 

Proof of the lemma 12.11 

i) The recurrence relations (|2.21l) are evident from the expression (|2.20p in the 
definition 12.11 

ii) Given the natural number k, and according to the definition (|2.20[) with 

M = 2 

B^Oil./fc) =X-( m )X-(/i 2 ) + (fc+l)0i 2 A-(Ml) + fc^l^ _ (M2)+fc(fc+l^VlM2- 

(A.l) 

Using (|2.7p . it is straightforward to verify that B^ipi, H2) — B^ (m2> Mi)- 
Moreover, it is clear from (|2.21j) that 

B^ifii, ■ • ■ )M»i ■ • • ) Mm) = ^i-i(Mi) • ■ • , ^i-i)B^ +l 1 (fii, Mi+i)-BM_i_i'' (m»+2> • • •! Mm) 

= b m (Mi, • • • > Mt+i, Mt> • • • ) Mm), 
therefore (mi, ■ • • , Mm) is a symmetric function. 

iii) Given the natural number fc, a direct consequence of the expression (|A.1|) is 
that 

S«( Ml)M2 ) = B^Oii.w) +?( M2 i?f ) (Mi) +MiBi fc) (M2)). 

If the expression iii) of the lemma HOI is assumed to be true for some M > 2, 
we have, using (|2.2ip . 

< +1 (m)= (i?r i) (M (M+1) ) + eE^<-i(M (2 ' M+1) )) (^ M+fc " 1) (MM +1 )+eMM +1 ) 

a/ 

= BK(M) + CMM+i X) (M (M+1) ) (<(M W ) + CMM+i B£'_ 1 (/* (i ' M+1) )) 

8=1 

M M+l 

= < 1 1 ) (M)+^^M(' t(i) )+^ +1 4 ) (f (M+1) ) = <!V)+fE M,<(M W ). 

t=l i=l 

Proof of the formula ([2~2"5f of the lemma [231 

Let k be a natural number. From the relations in (|2.7I) and the definition 12. II it 
is immediately clear that 

^(ALfM^ ^'^^'W (A.2) 
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Assume that the expression (|2.23[) of the lemma 12.21 is true for some M > 1 
then 



[MA),< + »] = [ma),<V m+1) )] ^ M+fe) (MM + i) + <(M (M+1) )[MA),^ M+fc) (^M + i)] 

= £ ( 2 <+1 (A U f ! ^ <+1 + ^ + 1} ) (^M(M^^) - 2*)) 

M 

= 2 ^ < +1 (AU^ )) + ^ g (^D (M (0) f^fo^"-*)) + — 2k) 



= £ f ^- Mi + ) (M W ) (m^m+i(^;M (i) ) - 2*) , 

where the equation (|A.2I) and the relation (|2.2ip of the lemma 12.11 are used 
when appropriate. 

The proofs of the formulas (|2.24l) and (12.25)) are analogous to the demon- 
stration above. 
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